Symmetric and Asymmetric Coalescence of Drops on a Substrate 
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The coalescence of viscous drops on a substrate is studied experimentally and theoretically. We 
consider cases where the drops can have different contact angles, leading to a very asymmetric 
coalescence process. Side view experiments reveal that the "bridge" connecting the drops evolves 
with self-similar dynamics, providing a new perspective on the coalescence of sessile drops. We 
show that the universal shape of the bridge is accurately described by similarity solutions of the 
one-dimensional lubrication equation. Our theory predicts a bridge that grows linearly in time and 
stresses the strong dependence on the contact angles. Without any adjustable parameters, we find 
quantitative agreement with all experimental observations. 
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The eoaleseenee or breakup of liquid drops is a fun- 
damental process relevant for the formation of raindrops 
or sprays, Inkjet printing, or stability of foams and emul- 
sions [l|-|3[ . The initial stages of coalescence of two spher- 
ical drops has been characterized in great detail |4l4lCll|. 
After contact, a small liquid bridge connects the two 
drops and the bridge grows rapidly with time. Depend- 
ing on the viscosity of the liquid, the radius of the bridge 
grows as r ^ t (high viscosity) JSnSl , or r ~ t^l"^ (low 
viscosity, inertia dominated) [THQI, with a crossover de- 
pending on fluid properties and drop size [13] ■ 

In many cases, however, the coalescing drops are not 
freely suspended but are in contact with a substrate. 
Much less is known about the coalescence of such ses- 
sile drops. When looking from a top view (perpendicular 
to the substrate), the coalescence of drops on a substrate 
looks very similar to the case for spherical drops Q; yet 
the bridge dynamics is fundamentally different. Mea- 
surements for very viscous drops give a growth r ~ i^/^ 



11 , [13] , and even smaller exponents have been suggested 
13|. The challenge lies in the complications introduced 
by the presence of the substrate. First, the geometry 
of the drop is no longer a sphere with an axisymmctric 
bridge, but a spherical cap with a contact angle 9. As a 
consequence, a top view of the coalescence process is very 
different from a side view. Second, the wall slows down 



the liquid transport towards the bridge [11| and gives rise 
to the motion of a contact line 1J|. At present, it is not 



clear whether or not this contact line motion affects the 
initial stages of coalescence, and different predictions for 
the dependence have been reported [ll|-[l3|. 

In this Letter wc resolve the coalescence of viscous 
drops on a substrate by performing side view experi- 
ments, imaging parallel to the substrate (Fig. [1]). Our 
central finding is that the initial stages evolve by a self- 
similar shape of the bridge, with a linear growth of the 
bridge height /iq ^ t. The influence of the contact angle 
is studied in detail by considering drops with identical 
or different contact angles, resulting into symmetric or 
asymmetric coalescence [Fig. [ijbc)]. Theoretically, we 




FIG. 1: (a) Schematic of two coalescing viscous drops on a 
substrate, viewed from the side. The minimum height /io(i) 
characterizes the bridge size. The left-right contact angles 
Ol and Or can be different at the moment of contact. The 
horizontal displacement 2;o results from the asymmetry in the 
contact angles. (b,c) Typical frames of the experiments are 
shown for asymmetric contact angles (b) and symmetric con- 
tact angles (c). 



show that all experiments can be described quantitatively 
by a similarity solution of the lubrication equation. Re- 
markably, this one-dimensional approach quantitatively 
predicts the shape and time evolution of the bridge with- 
out adjustable parameters. Our results reveal that the 
rate of vertical growth scales with the contact angle as 
^ 9'^, the horizontal speed ~ 9^, and provide a new per- 
spective on previous top view measurements. 

Experimental setup. — The side view images of co- 
alescing drops in Fig. dbc) arc obtained by a digital 
video camera (Photron APX-RS) equipped with a mi- 
croscopic lens (Navitar 12x zoom lens), resulting in a 
resolution of 2 ^m/pixel. The camera recorded 12.5 
frames per second. The substrate consists of a hori- 
zontal microscope glass slide (Menzel pre cleaned mi- 
croscope slide, average roughness w 10 nm). The glass 
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FIG. 2: (color online) Symmetric coalescence, (a) Height 
of the bridge ho as a function of time after contact t, for 
drops witli 6l ~ 6r ~ 22° (77 = 12.2 Pa ■ s). Experiments are 
shown in red (•), the solid line is the prediction by Eqs. (|2I3|I . 
The dashed line represents the lower limit for spatial resolu- 
tion, (b) Rescaled experimental profiles at different times, 
H = h{x,t)/ho{t) versus ^ = x6/ho(t). The collapse re- 
veals self-similar dynamics at the early stage of coalescence, 
in agreement with the similarity solution obtained from our 
analysis (solid line). 



a range of contact angles 9l and Or between 10° to 67° 
at the time of coalescence. We consider both symmetric 
coalescence [9l = Oji, Fig. [TJc)] and asymmetric coales- 
cence [Ol 7^ dft,, Fig.[IJb)]. The spreading determines the 
initial conditions, but in all cases the spreading speed is 
much smaller than the bridge growing speed. Contact 
time is determined when there is a visual change, which 
happens before the bridge is thick enough to provide a 
reliable measurement. The dashed line in Fig.[2ja) shows 
this spatial resolution limit. 

Self-similar dynamics. — The dynamics of coales- 
cence is characterized by the growth of the bridge con- 
necting the two drops. Figure Ufa) presents the mini- 
mum height of the bridge, h^, ss, & function of time for 
a symmetric coalescence experiment {Ol = Or = 22°). 
At early times, we observe a linear increase of the bridge 
height, i.e. ho ^ t, while at later times the coalescence 
slows down. In these final stages the height of the bridge 
becomes comparable to the total drop size, which is typ- 
ically ^ 1 mm for all experiments. The very early stage, 
however, exhibits self-similar dynamics that is governed 
by a single length scale. This is revealed in Fig. [2jb) 
where the meniscus profiles, h{x,t), and the horizontal 
coordinate, x, are rescaled by ho{t). The scaled profiles 
at different times collapse onto a universal curve: the 
early stages of coalescence are characterized by a self- 
similar meniscus profile. The size of the bridge is simply 
ho, both in horizontal and vertical direction. The solid 
line is the theoretical similarity profile that will be de- 
rived below. 

Our experiments suggest that coalescence of drops on 
a substrate is governed by a similarity solution of the 
flow. To simplify the three-dimensional geometry of the 
coalescence, we assume that the flow is predominantly 
oriented from the centers of the drops towards the co- 
alescing bridge, as suggested by Ristenpart et al. [ll[. 
We therefore attempt a similarity solution based on the 
one-dimensional lubrication theory for viscous flows [isj : 



slide was further cleaned using ethanol and acetone, then 
submerged in an ultrasonic bath and dried with filtered 
nitrogen gas. The coalescing drops were made from sil- 
icon oils (Basildon Chemical Company Limited), with 
viscosity ?/ = 0.974 Pa ■ s or 12.2 Pa • s. which both have 
a surface tension 7 = 21- lO^'' N • m^^ and density 
p = 975 kg-m~'^. The silicon oils perfectly wet the 
cleaned glass slide {9eq ~ 0). 

The coalescence of two drops is controlled as follows. 
A first drop is deposited from the syringe on the sub- 
strate. Although the silicon oil perfectly wets the glass, 
the spreading of these high viscosity drops is very slow, 
with the liquid contact angle decreasing slowly in time. 
Subsequently, the glass plate is displaced by a manual 
translation stage and a second drop is placed next to the 
first one. By controlling the expelled volume of silicon oil 
and the time between the deposition of drops we achieve 
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Here, h{x, t) is the meniscus profile viewed from the side, 
T] is the liquid viscosity and 7 denotes the surface tension. 
This lubrication equation is valid for small contact angles 
and represents mass conservation: the second term is the 
surface tension-driven fiux of liquid towards the bridge, 
causing a growth of the bridge {dh/dt > 0). 

Consistent with our experiments, Eq. ([1]) has a simi- 
larity solution that imposes a linear time-dependence. 



9x 



h(x,t) =vfH(0, with — 

vt 



(2) 



where 'H{^) is the similarity profile of the meniscus 
bridge. Here we incorporated the contact angle 9 in the 
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scaling of x, such that the condition dh/dx = 9 trans- 
lates to Ti' ~ 1. The correct scaling of the coalescence 
velocity with 9 then turns out to be 



V 



j9^ 



(3) 



where y is a numerical constant that still needs to be de- 
termined. In combination with ([1]) and ([2]), this provides 
an ODE for the similarity profile H(^): 



V 



(4) 



In order to solve Eq. which is a fourth order ODE 
with one unknown parameter V, five boundary condi- 
tions are required. At the center of the symmetric bridge 



H{o) = i, n'io) = n"'{o) = 0, 



(5) 



while far away the profile has to match a linear slope of 
contact angle 9. For the similarity variables this becomes 



■H"(oo) = 0, 7^'(oo) = l. 



(6) 



The boundary value problem ((31111) uniquely determines 
the similarity solution for symmetric drop coalescence. 
It was solved numerically using a shooting algorithm, 
from which we obtained both the dimensionless veloc- 
ity, V = 0.818809, and the similarity profile As 
the infiuence of the contact angle was scaled out, the 
solution describes the coalescence for all contact angles, 
within the lubrication assumption of small 9. 

The similarity solution indeed provides an accurate de- 
scription of the coalescence experiments. The solid line in 
Fig. ^a.) is the prediction ((Sj without adjustable param- 
eters. The solid line in Fig. ^h) is the similarity profile 
?^(^) obtained from our analysis. The agreement between 
experiment and theory shows that the initial stages of co- 
alescence are accurately described by a one-dimensional 
lubrication model. As expected, the similarity solution 
breaks down at later times when the size of the meniscus 
bridge becomes comparable to the size of the drops. 

Asymmetric coalescence. — We further extend the 
theory to asymmetric coalescence, for which the contact 
angles 9l ^ 9r (FiglH). Without loss of generality, we 
assume that 9^ > 9fi, and scale the coordinates using 9^. 
Interestingly, the lack of symmetry induces a horizontal 
displacement of the meniscus bridge during the coales- 
cence process: the minimum of the bridge, xq, is pulled 
towards the lower contact angle {9u). This effect can 
be captured using a similarity variable that is co-moving 
with the bridge, of the form 



^ {x — ut) 



with 
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(7) 
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FIG. 3: (color online) Asymmetric coalescence, (a) Horizon- 
tal and vertical position of the meniscus bridge, xo{t) and 
ho{t), for asymmetric drops {6l ~ 46°, 6r = 13°, viscosity 
r] = 12.2 Pa ■ s). Blue (•) and red (■) markers are experimen- 
tal data for 2:0 and ho respectively. The solid and dashed lines 
are the predictions from the similarity solutions, (b) Rescaled 
experimental profiles at difi'erent times, H = h{x, t)/ho{t) ver- 
sus 5 — xaOL/hoit). The collapse reveals self-similar dynam- 
ics at the early stage of coalescence. The solid line is the 
similarity solution predicted by our analysis. 



The horizontal velocity of coalescence u scales with 9\, 
where [/ is a numerical constant. The vertical velocity 
still follows © with 9 ^ 9l- Inserting © in ((T]) yields 



n 



1 {n^n'")' 



0. 



(8) 



This fourth order ODE for 'H(^) now contains two un- 
known parameters, U and V , hence the solution requires 
six boundary conditions. The minimum of the bridge is 
still defined by %{{)) = 1, W {Q) = 1, but the symmetry 
condition on Ti,'" no longer applies. Instead, one has to 
impose 00) = H'^oo) = 0, while the asymmetric 

contact angles give 



7^'(_oo) = -l, -^'(00) 



i/Ol 



(9) 
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FIG. 4: (color online) Contact angle dependence of coalescence velocity, (a) Dimensionless vertical speed, V = ivq / {^9\) , as a 
function of Ou/Ol- (b) Dimensionless horizontal speed, U = iurj/{'y9\), as a function of 9r/6l- The horizontal speed vanishes 
for the symmetric case 6r/9l = 1, and displays a maximum around 6l/Qr ~ 0.5. Symbols correspond to 75 experiments, solid 
lines are predictions obtained from the similarity solutions. 



With these conditions the boundary value problem has 
a unique solution for each ratio of contact angles 9r/9l- 
This means that U , V and H(^) can again be determined 
numerically for arbitrary 9r/0l. 

Figure [3] compares theory and experiment for an asym- 
metric coalescence [Or/Ol — 0.25). The horizontal posi- 
tion of the bridge xo (blue circles) and the vertical posi- 
tion of the bridge (red squares) arc shown in Fig.O^a). 
These again evolve linearly in time with a well-defined 
velocity. The solid and dashed lines are the predictions 
pi7p , with prefactors U and V determined from the simi- 
larity solution. FigurelSl^b) confirms that the asymmetric 
experimental profiles indeed display self-similarity (sym- 
bols), in excellent agreement with theory (solid line). 

We finally consider the influence of the contact angle 
on the coalescence speed. Our theory suggests a uni- 
versal behavior when making the horizontal and verti- 
cal velocities dimensionless, according ioU — 3mj/{'y9^) 
and V = ivrj / [-^9^) . The results of 75 experiments are 
summarized in Fig. 2] and compared to the theoretical 
prediction. Indeed, we observe a good collapse of the 
data. An interesting feature is that the theory predicts 
an optimal horizontal speed around 9^/9^ ~ 0.5, which 
is verified experimentally (although the experimental ve- 
locities arc slightly smaller than expected). This maxi- 
mum horizontal velocity can be explained as follows. The 
asymmetry induces a bias in the pulling force of surface 
tension, which is more efficient for the smaller contact 
angle 9^,. However, the "lubrication effect" inhibits liq- 
uid transport when — > 0, as the viscous friction in the 
liquid increases for smaller angles. The combination of 
these two effects gives rise to an optimum ratio 9b,/9l. 

Discussion. — Our results unply that the initial co- 
alescence of drops on a substrate, which is manifests 
three-dimensionally, is described quantitatively by a one- 



dimensional model. This can be explained from the cross- 
section of the bridge perpendicular to our viewpoint, r, 
which is much larger than h^. Elementary geometry sug- 
gests r {Rha/9y/^ 1^, R being the footprint radius 
of the drop on the substrate, such that indeed r ^ ho 
at early times. Combining this geometric relation with 
([3]) gives a growth law for the top view size of the bridge, 
r ~ 0^/^(i?7t/77)^/^. The dependence ^ 9^/^ differs from 
previous predictions [12[, but agrees with the data col- 
lapse proposed by Ristenpart et al. It would be 
interesting to see whether the one-dimensional approach 
also applies for coalescence of low- viscosity drops, which 
are dominated by inertia rather than viscosity [l^ . 

Our findings also highlight the key importance of 
asymmetry on the coalescence dynamics. While here the 
asymmetry is due to the contact angles, a similar effect 
was found for merging drops with unequal surface ten- 
sions (such as water and alcohol), for which the coales- 
cence is strongly delayed by Marangoni forces [l7| . This 



will have a strong bearing on applications as Inkjet print- 
ing, for which such asymmetries arc encountered natu- 
rally due to spreading and evaporation of ink drops. 
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